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IMAGING LARGE SCALE STRUCTURE IN THE X-RAY SKY 
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ABSTRACT 

We present the first results from a wide solid angle, moderately deep Chandra survey of the Lockman 
Hole North- West region. Our 9 ACIS-I fields cover an effective solid angle of 0.4 deg^ and reach a depth 
of 3 X IQ-i^ ergcm-2 s"^ in the 0.4-2 keV band and 3 x 10^^^ ergcm^^ g-i ^he 2-8 keV band. The 
best fit logN-logS for the entire field, the largest contiguous Chandra field yet observed, matches well onto 
that of the Chandra Deep Field North. We show that the full range of the 'cosmic variance' previously 
seen in different Chandra fields is reproduced in this small region of the sky. Counts-in-cells analysis 
shows that the hard band sources are more strongly correlated than the soft band sources. 

Subject headings: cosmology: observations — large-scale structure of the universe — x-rays: diffuse 
background — galaxies: nuclei 



1. INTRODUCTION 

Recent Chandra and XMM observations have resolved 
over 85% of the 2-8 keV X-ray background (XRB) into 
discrete sources, presumably active galaxtic nuclei (AGNs; 
Mushotzky et al. 2000; Brandt et al. 2001; Tozzi et al. 
2001; Campana et al. 2001; Cowie et al. 2002; Giacconi et 
al. 2002; Hasinger et al. 2001). However, the nature of the 
sources is unclear, with many of the faint sources showing 
little or no optical activity (Barger et al. 2001, 2002; Horn- 
schemeier et al. 2001; Rosati et al. 2002). There is a large 
field-to-field variance in the source counts in the Chandra 
deep surveys in the 2-8 keV band (Cowie et al. 2002). This 
cosmic variance is generally believed to arise from the un- 
derlying large-scale structure (LSS) that is traced, in some 
fashion, by the Chandra sources. Clustering of the XRB 
sources in redshift space has been seen in both Chandra 
Deep Fields (Barger et al. 2002; Hasinger 2002), indicat- 
ing that LSS does exist in the XRB source distribution. 
From a cosmological point of view, AGNs are expected to 
be highly biased tracers of cosmic structure formation at 
medium to large redshifts. The XRB sources, which have 

10 times the areal density of optically selected AGNs, 
provide for the first time a sufficiently high density to be 
used as tracers of the LSS. 

The clustering at various scales of X-ray selected AGNs 
has been pursued previously using observations from 
HEAO-1 ( e.g. Barcons & Fabian 1988; Mushotzky & 
Jahoda 1992) and ROSAT (e.g. Soltan et al. 1997; Vihk- 
linin & Forman 1995). The 05*^ T" results based on deep 
pointings (Carrera et al. 1998) show a much smaller am- 
plitude of correlation than predicted from optical sam- 
ples, while an analysis of the i? 05^ T North Ecliptic Pole 
sample, which is shallower but covers a wider solid an- 
gle (Mulhs 2002), comes to the opposite conclusion. It 
has been claimed that at a flux level below 1.5 x 10~^^ 
ergcm~^s~^ (0.5-2 keV) the XRB exhibits no clustering 



(Soltan & Hasinger 1994). This contradicts the observed 
cosmic variance in the ultradeep surveys by Chandra and 
XMM. How to combine these seemingly disparate results 
is one of the new mysteries of the XRB. 

To allow a direct measurement of the auto-correlation 
functions (ACFs) of the XRB sources, and to optically 
identify these objects to obtain their redhshifts, one needs 
arcsecond spatial resolution, a wide contiguous field of 
view, and sufficient sensitivity. To achieve these goals 
we have performed a large solid angle, moderately deep 
Chandra survey of the Lockman Hole North- West region. 
This survey currently has a contiguous sky coverage of 
~ 0.4 deg^ and is sensitive to X-ray flux levels of 3 x 10^^^ 
ergcm-^s-i (0.4-2 keV) and 3 x 10"^^ ergcm-^s'^ (2- 
8 keV). Most of the XRB is resolved at these flux levels 
(Cowie et al. 2002). 

2. OBSERVATIONS AND DATA REDUCTION 

The survey covers the Lockman Hole North- West re- 
gion centered at a = 10''34", S = 57°40' (J2000). The 
region has very low Galactic column density {Nh = 
5.72 X lO^^cm"^ , Dickey & Lockman 1990) and was cov- 
ered by the deepest 170/zm ISOPHOT field observed from 
ISO (hereafter ISONW). The field has also been observed 
at 850/im with the SCUBA camera on the James Clerk 
Maxwell Telescope and at cm wavelengths with the VLA, 
as well as at optical and near-IR wavelengths with the 
Subaru (using the unique Suprime-Cam instrument) and 
Keck telescopes. The ISONW region will be intensively 
observed with SIRTF as part of the LEGACY program. 
The multiwavelength analysis of this field will be presented 
in subsequent papers. 

The observation consists of 9 ACIS-I pointings (labeled 
as IS0NW1-IS0NW9) separated from each other by - 10' 
to allow close to uniform sky coverage. ISONWl has an 
exposure time of ~ 70 ks while the other 8 pointings have 
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Fig. 1. — The cumulative number counts in the 9 fields (a), (b) and the whole field (c), (d). For each (a) and (b), the insert shows the 
number counts of each ISONW field at flux level of 2 X 10~^^ ergcm" ^s" ^ (soft) and 9 X 10~^^ ergcm" -^s" ^ (hard). The solid line is the 
logN-logS of CDF-N. Field numbers 10 and 11 represents the CDF-N and CDF-S value respectively. The highest counts in the 9 fields are 
about 2.3 (soft) and 3.3 (hard) times higher than that from CDF-S, which has the lowest normalization among all the deep surveys (Cowie et 
al. 2002). For the hard band, number counts of 8 of the 9 fields are ~ Icr about the mean, but ISONW5 is much lower. For each (c) and (d), 
cumulative number counts for the whole field (diamonds with solid lines representing the Poisson error) is compared with that from CDF-N 
(solid histogram), CDF-S (dash-dotted histogram) and SSA13 (dashed line). The inserts shows the sensitive sky area vs flux. 



exposure times of ^ 40 ks. 

The data wore reduced with CI AO 2.2.1 and CALDB 
2.15. The CXC Science Threads were followed in data 
preparation. The resulting event lists were binned into 2 
energy bands, the soft (0.4-2 keV) and the hard (2-8 keV). 

Point sources were detected for each pointing in both en- 
ergy bands with wavdetect within the CIAO package. The 
wavelet scales of square root series 1, v^, 2, 2v^, 4, A^/2, 8 
and a false detection threshold of 1 x 10~^ were used. 
Spectrally weighted monochromatic exposure maps were 
created, assuming a power law with photon index of 1.2 
for the hard band and 1.4 for the soft band (Mushotzky 
et al. 2000; Barger et al. 2001). Count rates were 
converted to flux assuming the above power law spec- 
tra with only Galactic absorption. The conversion fac- 
tors are 4.74 x 10"^^ erg cm"^ (soft) and 2.34 x 10~" 
ergcm"^ (hard). The degradation of the ACIS low energy 
quantum efficency during the flight was corrected using 
the measurements of the ACIS team. 

The catalogs for all observations were merged. Source 
properties for objects detected in more than one observa- 
tion (due to the overlapping of fields) were taken from the 
field in which the source had the smallest off-axis angle. In 
the soft band, 431 sources were detected, and in the hard 
band, 278. The combined catalog contains 554 sources. 

3. ANALYSIS AND RESULTS 



Since the cfFcctive area decreases and the point spread 
function (PSF) increases with off-axis angle, the sensitiv- 
ity is not uniform across the field. We quantified this 
with Monte-Carlo simulations. First we constructed back- 
ground maps by removing the wavelet detected sources 
from the observed images and filled the holes with Pois- 
son noise sampled from regions surrounding the sources. 
Sources with fluxes drawn from the LogN-LogS derived 
from the Chandra deep flelds (Cowie et al. 2002; Garmire 
2002) wore generated and distributed uniformly within 8' 
from the aim point. PSF images of each source were made 
using mkpsf. The exposure maps described in the pre- 
vious section were then applied to simulate the effect of 
vignetting. The simulated sources were added to the back- 
ground maps to create simulated images. About 100 sim- 
ulations were performed for each band and exposure (for 
details see Yang et al. 2003). Above flux thresholds of 10 
cts/exposure (soft) and 12 cts/exposure (hard), the detec- 
tion is complete and the derived fluxes are consistent with 
the input values. We use the flux unit in cts/exposure 
because we found the wavdetect exposure-corrected count 
rates threshold in such unit only depends weakly on the 
exposure time in our observations due to the low back- 
grounds. We obtain the same results by assuming a uni- 
form background and the wavdetect detection thresholds 
described in Freeman et al.(2002). Within 4.5' of the aim 
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point, the detection is complete above 5 cts/exposure in 
the soft band and 7 cts/exposure in the hard band. 

Above flux thresholds of 1.2 x 10~^^ ergcm"^ s~^ (soft) 
and 8 x 10~^^ ergcm~^s~^ (hard), the source catalog in 
the combined field within 8' from each aim point is com- 
plete. This defines a flux limited sample that contains 115 
hard band and 298 soft band sources. The solid angle 
covered by the complete sample field is 0.33 deg^. 

We constructed the cumulative number counts (LogN- 
LogS) for each of 9 fields, as well as the whole merged 
catalog (Fig. 1). The Eddington bias found at these flux 
levels from the simulations is small and can be ignored. 



DEC ( deq) 



630, S'o = 2 X 10 



-15 




10h36m 10h34m 10h32in 

RA ( J2 ) 

Fig. 2. — The density map of the flux Hmited samples (region out- 
lined with a dashed line). The maps have been adaptivcly smoothed 
using fadapt in FTOOLS with a 6 counts kernel counts threshold. 
The gray scale map represents the source density in the soft band. 
The contours shows the source density in the hard band. The 7 con- 
tour levels are proportional to the logarithm of the source density. 

Fig. la and lb show comparisons of the LogN-LogS in 
the 9 pointings. Although the overlapping of fields could 
minimize the differences between them, this is the simplest 
way to demonstrate the cosmic variance because exactly 
the same procedures were used in each field and 8 of the 9 
fields have virtually identical exposures. Variance is seen 
in both the soft and hard bands. The full range of vari- 
ance seen in the Chandra deep surveys published to date 
is reproduced in this set of contiguous fields. 

The LogN-LogS of the whole field is compared with 
that from Chandra deep field North (CDF-N; Brandt et 
al 2001), the Chandra Deep Field south (CDF-S, Rosati et 
al. 2001) and SSA13 (Mushotzky et al. 2000) in Fig. Ic 
and Id. Our LogN-LogS connects smoothly with that of 
CDF-N, which has the highest normalization of all the 
published Chandra deep fields. 

The LogN-LogS distribution of the combined fields is 
modeled with power laws in the form of iV(> S) = 
A[S / So)~°' , using the area weighted maximum likelihood 
method (Murdoch, Crawford & Jauncey 1973). For the 
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0.23. For the hard band above 1x10 erg cm ^ s 



0.72 ±0.18; above 1 x 10"" 
= 1 X 10-^^,anda = 1.58 ± 
. the 

359, 5o = 1 X 10-'^^, anda = 



best fit parameters are A 
1.56 ±0.16. 

To visualize the source distributions, we used adaptive 
smoothing to create density maps of the flux limited sam- 
ple (Fig. 2). Structure is visible in both the soft and hard 
band maps, but the hard band sources are more clustered. 
To test whether the observed over-density could arise from 
Poisson fluctuations, we employed the likelihood test de- 
scribed in Carrera et al.(1998). Comparing with 10000 
simulated samples which are Poisson distributed, we found 
97.67% of the soft and 99.99% of the hard band simulations 
had better likelihood than the observations. The signifi- 
cance of clustering is therefore 2cr (soft) and 4cr (hard). 

Using the source distribution in cells tiling the field 
(counts-in-cclls), we can estimate the correlation scales 
of the sources. The variance of counts-in-cells defined as 
/Z2 = {{N — NY), where N is the mean counts in the 
cell, is directly related to the angular correlation function 
(Peebles 1980) by 



112 = N + 
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where ft is the cell size. The first term is the Poisson 
fluctuation. For correlation functions with a power-law 
form w{9) = (6/6oY^^ (where 7 is the power law index 
of the spatial correlation function) and square cells with 
size = 9x0 deg^, the integration can be obtained as 
(Totsuji & Kihara 1969; Lahav & Saslaw, 1992), 



2 _ 



(2) 



where is the normalized is a function of 7 

and is calculated numerically. We calculate cr^ for square 
cells of various sizes that tile the whole fleld. By fitting 
the cr^ — relation we should be able to estimate 60 and 7. 
We found the present data cannot constrain both param- 
eters accurately. By fixing 7 = 1.8, the "universal" slope 
measured in galaxies and in groups and clusters of galaxies 
(Bahcafl 1988), and minimizing x^, we found ^0 = 40±11" 
and ^0 = 4 ± 2" for the hard and soft band sources respec- 
tively. While the soft band sources agree very well with the 
angular correlation scale previously seen in ROSAT sur- 
veys (Vikhlinin & Forman, 1995), the hard band sources 
arc much more strongly correlated. 

The striking difference in clustering between the soft and 
hard band sources indicates the hard sources which are 
not detected in the soft band are highly clustered. About 
60% of the hard-band-only sources lie in overdense regions 
which form a 'band' connecting the 'lumps' on the west- 
ern and eastern side of the field (Fig. 2). This band in- 
cludes only about 1/3 of solid angle of the whole field. The 
counts-in-cells analysis (Fig. 3) also indicates that these 
hard-band-only sources have larger correlation scales than 
the soft band or hard band sources. 

4. DISCUSSION 

ASCA observations have shown that the rms variance of 
the 2-10 keV XRB on a scale of 0.5 deg^ is ~ 6%(Kushino 



soft band, between 2 x 10 and 10 erg cm ^ s ^ , ^4 = et al. 2002). With a sky coverage of 0.4 deg^ and a depth 
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that resolves > 50% of the hard band XRB, the norniaUza- 
tion of LogN-LogS derived from our observations sliould 
be very close to the "true" value. The fact that the ob- 
served LogN-logS connects onto the CDF-N field at low 
fluxes then indicates that > 90% of the XRB is resolved 
using the ASCA/ROSAT XRB normalization (Chen et al. 
1997). The main uncertainty is the normalization of XRB 
itself. 

The LSS seen in our field has reproduced the cosmic 
variance observed previously in deep field surveys, ft is 
noticeable (Fig. lb) that on scales of a Chandra field the 
variance is demonstrated as holes rather than lumps, i.e., 
the number counts in 9 fields are close to the mean, while 
only 2 field have very low value (IS0NW5 & CDF-S). This 
indicates the existence of voids in the X-ray LSS, which 
should, within a factor of a few, be of the same angular 
size as an ACIS-I field. Most of the hard-band-only sources 
cluster in relatively small regions, which may be topologi- 
cally connected. If this is not a result of a projection effect, 
then it is the first time a wall-like structure has been seen 
in the X-ray sky. 

We have found more variance in the hard band than 
in the soft band, consistent with the lower mean redshift 
of the hard X-ray sources. However, the difference in the 
mean redshift cannot account for most of the large differ- 
ence in the soft and hard band angular correlation length 
found by the counts-in-cells statistic. The mean redshift 
found previously for the Chandra/XMM deep fields sources 
is ~ 0.8. Unless most of the hard band sources are at very 
low redshift, which does not seem likely because we see no 
low redshift spikes in the redshift distribution of deep field 
surveys, the relatively low redshifts for the hard X-ray se- 
lected sources could not account for the factor of 10 larger 
differences in correlation scales seen in these sources. The 
similarity of the redshift distributions for the sources found 
in CDF-N and in an observation with a similar exposure to 
ours (Castander et al. 2003), indicates that the luminos- 
ity functions have likely been sampled below the "knee" 
for all redshifts in our observations. Thus we feel that the 
stronger correlation function for the hard sources seen in 
our data is likely due to their stronger spatial correlation 
rather than a redshift effect. 

The large variance in the Chandra source counts indi- 



cates that they must be highly biased tracers of matter. 
Direct calculation of the bias will require redshifts and 
an understanding of how the luminosity function changes 
with redshift, but it is already clear that the Chandra 
sources show much more variance than galaxy counts (Co- 
hen et al 2000) at similar optical magnitudes. 
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Fig. 3. — Normalized variance cr^ measured with different sized 
square cells which tiles the whole field. Squares; the hard band 
sources; triangles: the soft band sources. Crosses; the hard band 
only sources. The errors are estimated via a boots-trap technique. 
The lines shows the best fit variance with 7 = 1.8 (fixed) and 
do = 40" ± 11 for the hard band, 6*0 = 4" ± 2 for the soft band 
sources. The hard-band-only points are not fitted due to the small 
number of sources. 
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